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We demonstrate that superconductors which break time-reversal symmetry can exhibit thermo¬ 
electric properties, which are entirely different from the Ginzburg mechanism. As an example, we 
show that in the s-\-is superconducting state there is a reversible contribution to thermally induced 
supercurrent, whose direction is not invariant under time-reversal operation. Moreover in contrast 
to Ginzburg’s mechanism it has a singular behavior near the time-reversal symmetry breaking phase 
transition. The effect can be used to confirm or rule out the s + is state, which is widely expected 
to be realized in pnictide compounds Bai_ cc K a ,Fe 2 As 2 and stoichiometric LiFeAs. 
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Thermoelectric effects in superconductors, and their 
fundamental importance, were first discussed 70 years 
ago by Ginzburg [1, 2]. Currently there is a revival of 
interest in this topic [3-6]. These effects originate in 
charge transfer by thermal quasiparticles [7], which is 
compensated by the counterflow of superconducting cur¬ 
rent j s = — 5 n VT (where h n is the normal component 
thermoelectric coefficient). As it is determined by the 
dissipative normal current, such a thermally induced su¬ 
percurrent is irreversible since j s changes sign under the 
time-reversal transformation while b n and VT remain 
invariant. We point out here, that multicomponent su¬ 
perconductors which break the time-reversal symmetry, 
have entirely different reversible thermoelectric response. 
Namely, there is a generic contribution to the thermally 
induced supercurrent, whose direction is not preserved 
under time-reversal operation. Besides, that contribu¬ 
tion can be orders of magnitude larger than the ordinary 
one and, in contrast to Ginzburg mechanism, be present 
even at low temperatures. More precisely the new con¬ 
tribution we describe is important in the vicinity of the 
time-reversal symmetry breaking phase transition, that 
can occur at temperatures much lower than T c , where the 
usual contribution due to thermal quasiparticles is typ¬ 
ically extinct. To illustrate this new effect we consider 
time-reversal symmetry breaking states which are widely 
expected to exist in iron pnictide compounds [8-11]. 

A fundamentally interesting state that can appear in 
multicomponent superconducting systems is the so called 
s + is state. In addition to the usual gauge symme¬ 
try 17(1), it is characterized by a broken time-reversal 
symmetry (BTRS). Not only is it the simplest BTRS 
extension of the most abundant s-wave state, but also 
it is expected to arise from various microscopic physics 
[8, 9, 12-15]. The s + is state can as well be fabricated 
on demand on the interfaces of superconducting bilay¬ 
ers [16]. Recently it was demonstrated that such physics 
very likely occurs in strongly hole doped Bai_ x K x Fe 2 As 2 
[8, 15] as well as in stoichiometric compound LiFeAs [9]. 

In iron pnictides, the s + is BTRS state originates from 


the multiband character of superconductivity and sev¬ 
eral competing pairing channels, as shown schematically 
on the inset in Fig. 1(a). A typical band structure of 
Bai_ ;E K ;r Fe 2 As 2 consists of two hole pockets at the T 
point and two electron pockets at (0,7r) and (7r,0). The 
superconducting coupling here is dominated by the inter¬ 
band repulsion between electron and hole Fermi surfaces, 
as well as between the two hole pockets at T. Such a sys¬ 
tem can be described by a minimal three-band model, 
assuming that the order parameter components are Ai 
for electron pockets and A 2,3 for the hole bands at T 
[8]. The competition of repulsion forces results in an in¬ 
trinsically complex order parameter whose components 
in each band Ap = \Ak\e t0k (where k = 1,2,3 is the 
band index) possess non-trivial (frustrated) ground-state 
phase differences Okj = — 0j which are neither 0 nor n 

[8, 12-15, 17]. The ground-state order parameter of the 
s-\-is phase is not invariant under the discrete (Z 2 ) time- 
reversal transformation associated with complex conju¬ 
gation T(Afc) = AThus it has a broken U( 1) x Z 2 
symmetry [17] and below, s =b is will stand for states 
that are related by time-reversal transformation T. One 
of the interband phase differences, e.g. 0 12 can be used 
as the order parameter that characterizes the Z 2 phase 
transition. 

The BTRS states with phase frustration ^ 0, tt 
imply the existence of persistent “intrinsic Josephson-like 
currents” between the three bands, as shown schemati¬ 
cally in Fig. l(b,c). These currents have opposite “di¬ 
rections” for the two ground states with opposite values 
of interband phase differences. Although the underlying 
microscopic physics is not exotic, nor is any fine tuning is 
required to form such a state, no experimental observa¬ 
tions of such BTRS states have been reported. The rea¬ 
son being the major challenges to distinguish the s + is 
state from its time-reversal invariant s± and s ++ cousins 
by conventional methods. Indeed, the probe of relative 
phases between components of the order parameter in 
different bands is a non-trivial task. Few proposals have 
recently emerged for the indirect observation of BTRS 


2 





Figure 1. (Color online) - (a) Schematic view of the band structure in hole-doped iron pnictide compound Bai_ x K cc Fe 2 As 2 . It 
consists of two hole pockets at F shown by open circles and two electron pockets at (0, 7 r) and (7r, 0) shown by ellipses. The s+is 
state is favored by superconducting coupling dominated by the interband repulsion between electron and hole Fermi surfaces, 
as well as between the two hole pockets. Panels (b,c) schematically show the two degenerate ground-states s ± is. The phases 
of the order parameter components A k — \/Xk\e l ° k are represented by vector diagrams. Circular arrows show the directions 
of interband Josephson current supported by non-trivial relative phases Okj ^ 0, 7r. (d) Sample whose ends are maintained 
at different temperatures The thermophase effect appears due to the temperature-dependent interband phase difference 

Okj = 0kj(T). The total thermally induced current j will have opposite directions in s =t is states. Panel (e) sketches a closed 
circuit in spatially inhomogeneous bulk superconductor with two branches characterized by different values of the thermophase 
coefficients 7 ^ 1,2 ^(T). The junctions between branches have different temperatures T\^ and the temperature bias generates 
magnetic flux through the circuit (7). 


signatures in collective mode spectrum [15, 17-22], un¬ 
usual topological defects [23, 24] or spontaneous currents 
near impurities in samples subjected to strain [25]. 

The challenges associated with the observation of the 
s + is state and its relevance to pnictides motivate our 
choice of this kind of material for demonstration of the 
unusual thermoelectric properties of BTRS superconduc¬ 
tors. Below we show that in such systems, there is a new 
mechanism for thermal generation of supercurrent. It 
is not related to normal thermal currents and leads to 
thermoelectric properties which strongly differ from that 
suggested by Ginzburg, in usual superconductors. Al¬ 
though the suggested mechanism does not eliminate the 
usual thermoelectric effect, it will be shown to dominate 
near the BTRS phase transition thus providing an exper¬ 
imental test for this unconventional state. 

The key idea behind the proposed thermoelectric effect 
is that, due to the generically temperature-dependent in¬ 
terband phase differences, i.e., Okj = @kj{T ), a tempera¬ 
ture bias generates phase gradients of condensate compo¬ 
nents. Assuming that temperature gradients are small, 
so that the order parameter is determined by the local 
temperature, yields the relation 

V0* = V(A)+ 7 fc(T)VT, ( 1 ) 


effect. In multiband superconductors, the total cur¬ 
rent has contributions from each band j = J2 k jk- In 
units where h = c = 1 , the partial currents read as 
jk = (V0 fc -27rA/$o)/8e7rA2, where X k are characteristic 
constants associated with contributions of a given band 
to the London penetration depth. A is the vector poten¬ 
tial of the magnetic field B = V x A and 4>o the flux 
quantum. Introducing the notation Qo = V(^o~27rA/ff>o 
and the London penetration length A 2 = 1/ A^ 2 , al¬ 

lows to write the total current as 


3 = 


Qo 

8e7rA| 


+E 

k>i 


A? - X 2 k 
24e^A|Af 




( 2 ) 


The first term here is a usual Meissner current while the 
second part is determined by the gradients of interband 
phase differences. It describes the charge transfer by 
counter-currents of several superconducting components. 
Observe that according to Eq. ( 1 ), the counterflow is gen¬ 
erated by the temperature gradient and contributes to 
the current according to the following generic expression: 


. = Qo+7(^)VT 
3 8e7rA|(T) 

7 (T) = A|(T)^ 7 ft (T)A7(T). 

k 


( 3 ) 

( 4 ) 


where ip 0 = Y?k=i °k/3, and 7 k (T) = lT,j^k d0 kj/dT 
are the thermophase coefficients. Superconducting cur¬ 
rents can be generated, as a results of this thermophase 


The essence of such a thermophase effect is schematically 
shown in Fig. 1(d), for the case of Qo = 0 so that, accord¬ 
ing to Eq. (1), the partial currents in each band are given 
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by jk = 7 /c(T)VT/ 8 e 7 rA^. Since the thermophase coeffi¬ 
cients are opposite for the s±is states: T"fk = — 7 /e, these 
thermally induced superconducting currents are sensitive 
to time-reversal symmetry transformation. That is, for a 
given temperature bias, the current directions are oppo¬ 
site in s + is and in s — is states. 

Deep in the bulk of a superconducting sample the to¬ 
tal current (3) should vanish: j = 0, which yields the 
condition for current compensation 

Qo = - 7 VT. (5) 

This compensation can be achieved locally even for in¬ 
homogeneous superconductor if the thermophase coeffi¬ 
cient has slow spatial dependence 7 = 7 (T, r). Inte¬ 
grating Eq. (5) along a closed path and assuming for 
simplicity that no vortices are trapped in the circuit 
(<f dl • V(/?o = 0 ), we get the thermally induced magnetic 
flux 

$ r =ydl-VT 7 (T,r). ( 6 ) 

Since the coefficient y(T, r) is in general spatially inho¬ 
mogeneous along the integration path, the expression ( 6 ) 
yields a finite value of 4 >t- 

Now consider a closed circuit with two branches that 
are made up of different superconductors, for example 
the Bai_ :c K ;E Fe 2 As 2 compound with different doping lev¬ 
els. The thermophase coefficient y(T) in this case has 
a step-wise discontinuity along the circuit, determined 
by the values 7 ^ 1 , 2 ^(T) at the different branches. This 
is shown schematically in Fig. 1(e). Assuming that the 
junctions between branches have different temperatures 
Ti ? 2 , Eq. ( 6 ) yields an induced magnetic flux through the 
circuit given by 

^ £* dT[ 7 « (T) - 7 < 2 > (T)}. (7) 

The maximal possible magnitude of can be estimated 
by considering one of the branches to be in the time- 
reversal invariant state, so that e.g. 7 ( 2 \T) = 0 and 
Ai < A 2 . Then, from Eq. (7), we obtain 0 ~ 

50i 2 /47t where SO 12 = # 12 (^ 2 ) - # 12 (?i)- 

The resulting thermally induced flux can have a gi¬ 
ant magnitude compared to that produced by a usual 
thermoelectric effect. Below, we will introduce a realis¬ 
tic microscopic model for the s + is state, to demonstrate 
that the interband phase difference can have a significant 
variation as a function of temperature max(c) 0 i 2 ) ^ 7 r. 
Therefore the resulting flux 4 >t ^ 4>o is much larger than 
the typical value of 4 >t ~ 10 -2< ho expected in usual su¬ 
perconducting thermoelectric circuits [1, 2, 26—28] 1 . 


The magnitude of quasiparticle thermoelectric effects in su¬ 


To provide a microscopic basis for this physics, we 
now proceed to calculate thermophase coefficients for the 
s-his superconducting state in Bai_ ;E K ;E Fe 2 As 2 . Within 
the minimal three-band model, it is parametrized by two 
pairing constants characterizing the strength of interband 
hole-hole Uhh and electron-hole Uhe repulsions (as shown 
schematically in Fig. 1). The self-consistency equation 
for the order parameter components A = (Ai,A 2 ,As) 
[ 8 , 15, 21] 

N d 

A = 2nTA y] F(u n ) , ( 8 ) 

n =0 

where F = (F 1 ,F 2 ,F 3 ) and F k (uj n ) = A k /^wl + |A fc | 2 , 
uj n = (2n+ 1)7 tT are fermionic Matsubara frequencies, T 
is the temperature. Here Nd = Qd/(2 ttT) is a cut-off at 
Debye frequency, and the coupling parameter matrix is 
chosen in the form 

( 0 Uhh Uhe \ 

Uhh 0 Uhe I 5 (9) 

Uhe Uhe 0 J 

which yields the phase frustration for certain parame¬ 
ters Uhh') u he • In Ibis case Ai ? 2 and A 3 describe the 
order parameter in hole bands and electronic pockets 
correspondingly. The order parameter has symmetric 
form A = e^°(|Ai|e 2012 , |A 2 |e —2,6,12 , |A 3 |), where (p Q is 
the common phase and the non-trivial interband phase 
difference 612 7 ^ 0, 7 r appears in the TRSB s + is state. 
The phase diagram, shown in Fig. 2, is obtained by solv¬ 
ing a nonlinear self-consistency problem (8,9). It demon¬ 
strates the line of Z 2 BTRS phase transition that oc¬ 
curs at the temperature T % 2 < T c . Within this model 
the thermophase coefficients in Eq. (3) are given by 
7 i(T) = — 7 2 (T) = —\dd\ 2 ldT and 73 = 0. Despite the 
symmetry between bands the net thermophase coefficient 
(4) is in general non-zero since the coefficients A& depend 
not only on the gap amplitudes but also on the kinetic 
coefficients. For example in the diffusive limit they are 
given by a standard text-book expression for the London 
length A ^ 2 = 47 r 2 cr/c | A/~| tanh(| A&I/ 2 T) where or are the 
normal state conductivities in each band which are differ¬ 
ent even if the superconducting pairing is symmetric be¬ 
tween hole bands. In the vicinity of T % 2 the thermophase 
coefficients diverge as y(T) ~ (T % 2 — T) -1 / 2 . In Fig. 2(b) 
the overall profile of the function ( T % 2 — T)y(T) is shown, 
where the prefactor ( T % 2 — T) is added to remove the di¬ 
vergence. The typical temperature dependencies of y(T) 
are shown in Fig. 2(c) for particular coupling parameters. 


perconductors is still a challenging problem, see discussion in 
Ref. [29]. It is beyond the scope of the present paper since we 
deal with unconventional thermoelectric effect, which is not re¬ 
lated to the flow of quasiparticles. 
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The square-root singularity dependence of thermophase 
coefficient (see Fig. 2c) is consistent with the fact that the 
effective potential for the phase difference gets soft near 
the Z 2 symmetry-restoring phase transition. However it 
is in strong contrast with the quasiparticle contribution 
which has a smooth thermally activated behavior [26, 27]. 

Hence the thermophase coefficients are large near T % 2 
and the discussed thermoelectric effect is generically the 
dominating one at low temperatures in the fully gapped 
BTRS s + is state. The magnitude of the total variation 
of the relative phase 0 12 as a function of temperature can 
be estimated from the typical dependence in Fig. 2(c). 
There 8 O 12 ~ 0.47T, so that the thermally generated flux 
is ~ 0.24>q which is in general much larger than the 
theoretically predicted value for quasiparticle contribu¬ 
tion. 

To draw a more convincing comparison let us calcu¬ 
late the thermoelectric fluxes in the BTRS state TRS 
and in the usual time-reversal symmetric superconduc¬ 
tor <$>Jp RS generated in the thermoelectric circuit accord¬ 
ing to Eq.(7). For simplicity we assume that one of the 
branches is passive 72 = 0 and considerer the tempera- 





Figure 2. (Color online) - (a) Phase diagram of the three- 
band superconductor described by the microscopic model (see 
text). Solid curves show the critical temperature T% 2 of the 
BTRS phase transition, (b) Color scale plot of the ther¬ 
mophase coefficient (Tz 2 — T)y(T). The prefactor (T% 2 — T ) is 
added to regularize the square-root singularity of y(T) at the 
BTRS phase transition, (c) The temperature dependencies 
of 7 (T) along the lines such as shown by arrow in the panel 
(b). From right to left Uhe/uhh = 1-05; 1.15; 1.25; 1.35. 
(d) The ratio of magnetic fluxes generated due to the uncon¬ 
ventional and conventional thermoelectric effects, & RTRS and 
correspondingly. S n is a dimensionless normal state 
Seebeck coefficient which can be of the order 1 but typically 
is much smaller. The other system parameters are explained 
in text. 


ture bias T 2 — T\ = 0.1T C which is of the order of 1 — 5 K 
in iron pnictides. To calculate a usual quasiparticle con¬ 
tribution we use the result of Refs. [26, 27]: 7 = 47 tA \r\ 
where r\ = 0.6 a n S n (T/T c ) f^ o/T dyy 2 cosh - 2 (y/ 2 ). Here 
Aq is the smallest gap in the system, cr n is the nor¬ 
mal state conductivity and S n = ( e/kB)S n where S n 
is a Seebeck coefficient just above T c , ks is Boltzmann 
constant and e is the electron charge. We use a diffu¬ 
sive model for A& parametrized by kinetic coefficients 
<73 — (J 2 = olg\ where the coefficient a ~ \ measures 
an asymmetry between hole pockets. Note that within 
our simplified model of equal coupling constants for the 
hole bands the overall magnitude of the unconventional 
thermoelectric effect is proportional to the asymmetry 
of kinetic coefficients — a 2 = (a — l)cri, which we 
assume below to be rather small a = 0.9. For such pa¬ 
rameters the typical temperature dependencies of the ra¬ 
tio S n &BTRs/&TRS are shown in Fig. 2(d). From this 
plot one can see that even for quite a large values of 
Seebeck coefficient S n ~ 100 fiV/K and hence S n ^ 1 
observed in some iron-oxipnictide compounds [30] we ob¬ 
tain that the unconventional thermoelectric signal is in 
fact much larger than the quasiparticle contribution es¬ 
pecially if Tz 2 T c , see black curve in Fig. 2(d). 

The usual thermoelectric effect can be significant at 
low temperatures if quasiparticle density remains finite, 
e.g., in gapless superconductors. However, we would like 
to emphasize that even in this case the suggested effect 
can be separated from quasiparticle contribution since it 
breaks the time-reversal symmetry. That is, the direction 
of thermally induced magnetic field is different in s + is 
and s — is states while the usual thermoelectric effect is 
time-reversal invariant. 

Let us now consider several characteristic examples 
which should arise in concrete experimental set-ups based 
on circuits containing BTRS s + is superconductors. To 
that end, we use the three-component Ginzburg-Landau 
(GL) theory (10), although the effects will be similar in 
the microscopic formalism. The GL expansion, derived 
from the microscopic three-band model (see Appendix 
for details), reads as 



+ EE %j|Afc||Aj| cos 9kj + — , 

k=lj>k 


Pk 

2 


|A ,| 4 


( 10 ) 


where j, k are band indices, and are diffusion con¬ 
stants. We model temperature dependence of the coeffi¬ 
cients as ak — (T/Tfc — 1 ) (a^ and being charac¬ 
teristic constants, further details on the derivation and 
choice of parameters are given in Appendix ). The nu¬ 
merical studies of the GL equations are performed within 
a finite-element formulation [31] and using a non-linear 
conjugate gradient algorithm, with the standard condi- 
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Figure 3. (Color online) - Thermally induced magnetic field 
on a domain-wall between s + is and s — is phases. Two 
faces of the stripe are maintained at different temperatures, 
thus resulting in a thermal gradient , here ST = 0.2T c , along 
the domain-wall. The s + is/s — is domain-wall generates 
the magnetic field B t- In the absence of the temperature 
gradient a domain-wall in s + is superconductor does not have 
any magnetic signature. 


which does not affect the total thermally induced mag¬ 
netic flux. For a homogeneous BTRS state in the absence 
of domain walls both normal and unconventional thermo¬ 
electric signals are absent. 

The non-trivial thermoelectric effect in s + is phases 
can be employed to obtain an advanced functionality 
in the practical applications of certain iron-based su¬ 
perconductors. For this purpose one can manipulate 
the symmetry of superconducting states by simultane¬ 
ously applying temperature bias and an external mag¬ 
netic field H ext which in general removes the degeneracy 
of s + is/s — is states. As a result it is possible to pre¬ 
pare a particular state by cooling the system through a 
BTRS transition. As an example we suggest a scheme 
of a memory cell based on the non-trivial thermoelec¬ 
tric properties of bimetallic rings shown in Fig. 4 with 
both branches made of s + is superconductor, but hav¬ 
ing different doping. The junctions between branches are 
maintained at different temperatures Tp 2 as sketched in 
Fig. 4(a). There, the s ± is states have the opposite val- 


tion that no current flows through the sample’s bound¬ 
aries. 

First, we consider thermoelectric response of a sample 
with a domain wall that separate s d= is phases, which 
are natural inhomogeneities expected to exist in realistic 
systems with spontaneously broken discrete symmetries. 
It is known that such domain walls do not carry mag¬ 
netic field at constant temperature [24]. We find that, in 
the presence of a temperature gradient along the domain 
wall, it generates a thermally induced magnetic field, as 
shown in Fig. 3. There, the thermophase coefficients 7 (T) 
have opposite signs ins+is/s — is domains. Therefore, in 
the vicinity of the interface between these, there should 
be net superconducting current and a thermally induced 
magnetic field Bt . Note that there is no phase winding 
in this flux-carrying configuration and the thermally in¬ 
duced flux B t is not quantized 2 . Since the condensate 
density near domain wall is non-homogeneous one cannot 
avoid additional signal produced by usual the thermo¬ 
electric effect according to Ginzburg’s mechanism. How¬ 
ever, it generates a dipole-like magnetic field distribution 


2 Magnetic flux penetration into a superconductor without phase 
winding/vortices is not restricted to three-band superconduc¬ 
tors. This can be seen from the general relation for magnetic 
flux obtained from the expression supercurrent <f> = j> Adi = 

§[ J /pi r ) + J2i=l..N l A i( r )| 2V6, i( r )/p( r )] (p( r ) bein § the total 

density). Even if there is no phase winding and even if we select 
the integration contour where J —> 0, the integral will in general 
be nonzero, in the presence of gradients of relative phases 6i — 0j 
and densities. 



Figure 4. (Color online) - A bimetallic ring consisting of 
two branches in the s =b is state, but having different chem¬ 
ical composition (see details in Appendix ). Both junctions 
are maintained at different temperatures Ti, 2 , so that there is 
a thermal gradient that varies linearly with the polar angle. 
This is sketched in panel (a). The resulting thermally induced 
magnetic field is shown on panel (c). In contrast to the usual 
thermoelectric effect its overall sign depends not only on the 
temperature bias but also on which of Z 2 BTRS ground states 
(i.e. s + is or s — is) is realized. This configuration supports 
nonzero flux 4 >t in the inner hole, which is more or less inde¬ 
pendent of the diameter of the inner hole. The total flux, on 
the other hand depends on the diameter of the sample (b). 
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ues of Bt 3 . In such a system the value of the Z 2 
index, which distinguishes s ± is states, can be consid¬ 
ered as a bit of information. Let us now describe the 
writing/reading protocols. First the external field H ext 
interacts with thermally induced currents and removes 
the s ± is degeneracy. Hence the particular s + is or 
s — is state can be selected when such a system is cooled 
through the BTRS phase transition, in the presence of 
both external field and temperature bias ST = T 2 — T\. 
After cooling, both the external field and temperature 
bias can be removed, leaving the system in the ground 
state characterized by the Z 2 index sgn (STH ext ). The 
read-out protocol, i.e. the measurement of Z 2 index, can 
be implemented by applying the temperature bias (in the 
absence of external field) and measuring the direction of 
induced magnetic flux. 

To conclude, we demonstrated that multicomponent 
superconductors with broken time-reversal symmetry, 
and in particular the s + is state, feature in the vicinity of 
BTRS transition, a giant thermoelectric effect of princi¬ 
pally different nature than that in single-component su¬ 
perconductors. It originates in thermally induced inter¬ 
component counter-flows, in contrast to the counter-flows 
of normal and superconducting currents in the mecha¬ 
nism originally discussed by Ginzburg. Although related 
effects should be present in various multicomponent su¬ 
perconductors, along with multicomponent superfluids, 
we focused on the s + is states where the effect is generic 
(irrespective of the model). 

The work was supported by the Swedish Research 
Council grants 642-2013-7837, 325-2009-7664. The com¬ 
putations were performed on resources provided by the 
Swedish National Infrastructure for Computing (SNIC) 
at National Supercomputer Center at Linkoping, Swe¬ 
den. 


3 The unconventional pattern of magnetic field induced by 
counter-flows is related to the general property of non-Meissner 
contribution to electrodynamics, in multicomponent supercon¬ 
ductors. The essence of the physics leading to the Meissner 
and, in a broader sense, to the Anderson-Higgs effect is that the 
magnetic field is described by a massive vector field theory. As 
shown in Refs. [32, 33], multicomponent systems cannot gener- 
ically be described by a massive vector field theory. Rather, it 
can be mapped onto a massive vector field theory coupled to 
another field associated with the charge-carrying counter-flows. 
The later terms lead to contributions which are typically negligi¬ 
ble or unobservable in ordinary multiband superconductors due 
for example to interband couplings. By contrast as discussed in 
this paper, the thermoelectric effect in the s + is state, makes 
the interband phase fluctuations propagate into the bulk. This 
produces the unconventional magnetic signatures. 


Ginzburg-Landau expansion for multiband 
superconductors 

To derive the Ginzburg-Landau expansion (10), we use 
the microscopic self-consistency equation 

N d 

A = 2ttTA^F(0, (11) 

n=0 

where F = (F 1 ,F 2 ,F 3 ) and F k (uj n ) = A k /^wl + |A fe | 2 . 
uj n = ( 2 n+l) 7 rT, where n G Z is the fermionic Matsubara 
frequency, T the temperature, and Nd = f]^/(27rT) is 
a cut-off at Debye frequency. In the diffusive case the 
anomalous functions in each band can be expanded as 
follows: 

Aj AjIAjP D.fV-ieA)^ 

• (12) 

where Dj is the diffusion coefficient. Then we get for the 
summation over Matsubara frequencies 


oo 


27 tT ^2 = Go + r, 

n=0 

(13a) 

oo 

= 0 . 1 /T 2 , 

n=0 

(13b) 

oo 

2ttT^1/(2 W 2 ) = 0.4/T, 

(13c) 


n=0 


where r = (1 — T/T c ). We normalize A by T c and 
Co = \ZnDo/8f~c, where D 0 is some arbitrary diffusion 
constant. 

(G 0 + r — A -1 ) A — 6| A| 2 ■ A = D ■ (II 2 A) (14) 

where b = 0.1, II = V — icA. A = (Ai, A_>. A 3 ), | A| 2 = 
(|A 1 | 2 ,|A 2 | 2 ,|A 3 | 2 ), and D=(D 1 ,D 2 ,D 3 ) where Dj = 
Dj/Do. The diffusion coefficients can be different. In the 
diffusive case the mixed gradient terms are absent. 

Let us consider the case of the intraband dominated 
pairing which can be described by a three-component GL 
theory in the vicinity of T c . This regime is described by 
the following coupling matrix: 

( A ~Ve \ 

A = ~Vh A -rje (15) 

V -Ve -rje A / 

where r/ e , rjh <C A. The critical temperature is determined 

by the equation Go = min(A^ 1 , A^" 1 , A 3 x ), where N[\ = 

( 2 A - 1 i h ± a/ 87? 2 + ^)/[ 2 (A 2 - A r) h - 2 r? 2 )] and Ag 1 = 
1 /(A + rjh) are the positive eigenvalues of the matrix 

( A 2 -?? 2 Ve+ X Vh 7 ie(^ + Vh)\ 

A -1 = X I rjl + Xrih A 2 - ry 2 r] e {X + %) (16) 

\Ve(X + Vh) Vei X + Vh) A 2 — r)\ J 
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where X = 1/[(A 2 — Xrjh — 2 ^) (A + 77 ^)]. Since we assume 
that Tj e ,h > 0 and rj e ^ <^C A the critical temperature is 
given by the smallest eigenvalue Go = 1 / (A + 77 ^) so that 

( ax ax a 2 \ 

ax ax a 2 1 (17) 

a 2 a 2 a 3 ) 

where 

a\ = {ril + \-q h ) / X (18a) 

a 2 = r) e {\ + rj h )/X (18b) 

«3 = (2 ??e ~Vh + *nh)/X . (18c) 

For example, for A = 1 and r\ e = 0.1, rjh = 0.2 we get 
ax = 0.2244, a 2 = 0.1282 and <23 = 0.1923. For A = 1 
and T] e = rjh = 0.1 we get ax = a 2 = as = 0.1136. 

Then we get the system of GL equations 

(r — ai) A] — aiA 2 — a 2 As — 6Ai|Ai| 2 + DiII 2 Ai = 0 

(r — ai)A 2 — aiAi — a 2 As — 6A2IA2I 2 + D 2 Tl 2 A 2 = 0 

(r — as)As — a 2 Ax — a 2 A 2 — 6A3IA3I 2 + Dsll 2 As = 0 

where b = 0.1, II = V — ieA and Dj = Dj/Dq. The 

diffusion coefficients can be different. In the diffusive case 
the mixed gradient terms are absent. The free energy 
functional whose variations gives the Ginzburg-Landau 
equations reads as 

J 7 = |nA/e| 2 + ak\A k \ 2 + ^|A/e| 4 | 

k=1 
3 3 

+ E E Vk^AlAj + A*A k ] , (19) 

k= 1 j=k+ 1 

where /? fe = 6, 7712 = ai, hi 3 = V 23 = a 2 , a k = a° k (T/T k - 
1), and a° k = l- a k , T k = T c { 1 - a k ). 

Details of simulations 

In order to investigate the response to the thermal gra¬ 
dients, we minimize numerically the free energy (19), in 
zero external field, with the standard condition that no 
current flows through the sample’s boundaries. In sim¬ 
ulations we use the microscopically derived coefficients 
(18). The theory is discretized within a finite element 
formulation [31] and the minimization is performed us¬ 
ing a non-linear conjugate gradient algorithm. 

Details for the bimetallic ring 

In Fig. 4 of the main text we consider the case of a 
bimetallic that consists of two branches in the s + is 
state, but with different chemical composition. The dif¬ 
ference in the chemical composition is modelled by modu¬ 
lating the Ginzburg-Landau coefficients (18) in the lower 


branch to be 85% of their value in the upper branch: 
af °wer) = 0.85a^ upper) . Thus the cq vary stepwise while 
passing from one branch to the other. The interfaces be¬ 
tween both branches are maintained at different temper¬ 
atures T\ and T 2 . Here this is the whole interface that is 
maintained at a given temperature, and thus the temper¬ 
ature is a linear function of the polar angle. This set-up 
determines the spatially inhomogeneous coefficients ak 
and rjkj used in Fig. 4. 
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